Quantum Inductance and High Frequency Oscillators in Graphene Nanoribbons 
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Here we investigate high frequency AC transport through narrow graphene nanoribbons with 
topgate potentials that form a localized quantum dot. We show that as a consequence of the finite 
dwell time of an electron inside the quantum dot (QD), the QD behaves like a classical inductor at 
sufficiently high frequencies uj >50 GHz. When the geometric capacitance of the topgate and the 
quantum capacitance of the nanoribbon are accounted for, the admittance of the device behaves like 
a classical serial RLC circuit with resonant frequencies u ~ fOO — 900 GHz and Q-factors greater 
than 10^. These results indicate that graphene nanoribbons can serve as all-electronic ultra- high 
frequency oscillators and filters thereby extending the reach of high frequency electronics into new 
domains. 



I. INTRODUCTION 

The recent isolation of graphene fT'| , a two dimensional 
atomically thin crystal comprised of sp^ hybridized car- 
bon atoms, sparked an unprecedented amount of research 
activity aimed at understanding and exploiting the un- 
usual properties of this material. Early research efforts 
centered around understanding graphene's fundamental 
properties, and their applications for electronic devices. 
Graphene was shown to exhibit anomalous 0-01 and frac- 
tional [sl, [1| quantum hall effects, 7r-Berry phase and 
be capable of ballistic 043; and coherent [l^ trans- 
port. Furthermore, its ultra high mobilities, and the 
promise of realizing ballistic transport at elevated tem- 
peratures attracted the attention of device physicists and 
engineers, who soon showed that properties such as high 
room temperature mobilities [TTl, excellent thermal con- 
ductivity ^12i] . and unusually high mechanical durability 
[l3l |. render graphene to be the ideal material for single 
molecule gas detectors [ij, [l^, high density capacitors 
[l^ . and most notably, ultra-high frequency transistors, 
which were recently shown to be capable of 100 GHz 
operation [TtJ, and predicted to be capable of THz op- 
erating frequencies (l8l420( . Consequently, graphene is 
listed as one of the candidate materials for post-silicon 
electronics on the International Technology Roadmap for 
Semiconductors (2l| . 

Although DC and low frequency transport in graphene 
nanostructures is well understood j.22;], high frequency 
transport in these devices has not yet been thoroughly 
investigated. Furthermore, the sparse experimental mea- 
surements that characterize high frequency graphene de- 
vices are limited by the maximum frequencies which can 
be directly measured using commercial electronics. Con- 
sequently, the cut-off frequency in these devices can only 
be extrapolated from its 1 / /-gain plot [l^, H^l , but not 
directly measured. In this work, we examine AC trans- 
port in top-gated graphene nanoribbons (GNRs) using 
the Green-Kubo linear response formalism. We show 
that at sufficiently high frequencies, transport in GNRs 
is analogous to a classical RLC circuit, where the induc- 
tive component, which is quantum mechanical in origin, 
becomes dominant after a transition frequency ujtr- Fur- 
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Figure 1: (a) Schematic representation of the device and the 
biasing scheme, (b) equivalent classical circuit model, in which 
Cq is the quantum capacitance of the graphene nanoribbon 
and Lq is the quantum mechanical inductance of the top gate 
defined quantum dot. Cbg and Ctg are the geometrical ca- 
pacitances arising from the back and top gate coupling re- 
spectively. 



thermore, if this inductive behavior is coupled with the 
GNR's gate-tunable quantum capacitance, the resultant 
circuit can be utilized as an all electronic ultra high fre- 
quency oscillator. The ability to adjust quantum capac- 
itance of the CNR oscillator in-situ renders this device 
ideal for ultra high frequency all electronic switching and 
measurement applications. Furthermore, we propose a 
scheme in which GNR-based devices could be utilized to 
measure ultra high frequency signals, thereby overcoming 
current measurement limitations. 



II. THEORETICAL MODEL 

A schematic representation of the device under con- 
sideration and a biasing scheme are shown in Figure la. 
The device consists of a narrow graphene nanoribbon of 



width < 50nm, with Ohmic source and drain contacts 
and capacitively coupled back and top gates, which are 
electricaUy isolated from the GNR by two dielectrics with 
permettivities Si and £2- A thin metallic top gate is uti- 
lized to provide a locally tunable barrier via the applica- 
tion of an electrostatic bias Vtg, whereas a global bias 
may be applied to the backgate Vbgi which is used to 
modulate the Fermi level in the entire device. Further- 
more, the top gate potential is used to electrostatically 
define a quantum dot (QD). Similar devices are routinely 
fabricated (TgI . [25| , and transport in this structures has 
been studied in both DC f25l. l26l| and more recently low 
frequency AC |23| bias regimes. In our discussion we as- 
sume that only a single level of the QD is accessible to 
electrons in the source-drain bias window so that the dot 
can be characterized by a single energy level Eq and line 
width jL- 

In order to explore transport through this device in 
both DC and high frequency AC regimes, we utilize the 
Green-Kubo formalism to calculate the frequency depen- 
dent admittance Td {uj,x) through the QD using 

{lo,x) = j dyj d<e'(-+-)* (j(x,0,i(a:,0)) , 

-L/2 

(1) 

where L is the topgate length, j is the current operator, 
and e is a positive infinitesimal. The admittance for the 
top gate defined quantum dot, which we denote as F/j is 
evaluated in Appendix I. Furthermore, in order to couple 
the high frequency admittance of the QD to the GNR 
device, we utilize the method introduced by Wang et al 

2 9 9 9 

+ + (2) 

C^{lo) Co iVoit^^x) Cq' 

where, Co is the total geometric capacitance, and Cq 
is the quantum capacitance of the nanoribbon, which is 
proportional to the density of states at the Fermi en- 
ergy. The frequency dependent electrochemical capaci- 
tance C^(aj) is related to the admittance of the entire 
device by T{uj) = ~iLuCfj_{uj). 

A. Quantum Inductance 

In order to interpret the physical meaning of the com- 
plex QD admittance Td {uj,x) — > To (i^), we express it 
in terms of its real and imaginary components Td (i-^) = 
(w) + iSFu (u)) (see Appendix I for the complete 
expression). Similar expressions for the admittance are 
obtained in (27l . [28| using a nonequilibrium Green's func- 
tion approach. We further introduce a dimensionless pa- 
rameter C = 2 {Eq — fi) (see Appendix I) in order 
to characterize the width of the resonance, where fj, is 




-2-1 1 2 

Vtg (Volt) 




-4 -2 2 4 

Ftg (Volt) 



Figure 2: (a) Total gate tunable capacitance Ce of graphene 
using suspended top gate (black), 1 nm Si02 (blue), and 1 nm 
Hf02 (red) dielectrics; (b) a comparison between the quantum 
capacitance of bulk graphene (black) and that of a 2.5 nm 
graphene nanoribbon (purple) showing van Hove singularities 
on a suspended top-gate. 

the electrochemical potential of the source lead, and the 
difference Eq — defines the transport window created 
by the source-drain bias. In the limit a; 0^, we re- 
cover the DC Landauer conductivity 2e^ {h (l + C'^)) 
A plot of both real 5RF (w) and imaginary SF (w) com- 
ponents of the dynamic admittance are shown in Ap- 
pendix I. As can be seen in Fig. 7b, the sign of the 
imaginary component of the admittance becomes positive 
after a critical transition frequency wt/j, corresponding 
to a negative capacitance, or, equivalently, to inductive 
behavior. It should be noted that in quantum trans- 
port the sign convention is opposite to the one used in 
electrical engineering. Namely, the sign of the frequency 
component (e~*"*) is chosen such that 3F (w) > cor- 
responds to inductive behavior. Quantum inductance in 
various nanostructures has been previously investigated 
both theoretically [23,[2^,[33 and verified experimentally 

Phenomenologically, the inductive behavior of the QD 
can be understood by viewing the QD and topgate as 
forming a parallel plate capacitor. In classical parallel 
plate capacitors, charge accumulation on the device re- 
sults from the application of a voltage. However, even 
for a quantum coherent capacitor there exists an intrin- 
sic charge relaxation resistance Rq = /i/2e^ = 12.9 kil, 
in the limit of a single transport channel . Conse- 
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Figure 3: Gate and bias voltage dependence of tiie Fano Fac- 
tor of a GNR sk. uj = 800 MHz and (a) W/L = 25 and (b) 
W/L = 5, with L = 200nm. The se p arameters were chosen 
to match the experiments in Ref. (33. [34|. 



quently, the charge accumulation time is trc = {RqC)^^ . 
If the driving voltage is an AC signal, the charge accu- 
mulation will follow the voltage and will also reverse sign 
in due time . Ideally, in the absence of resistance, there is 
a 7r/2 phase lag between the current and the AC driving 
voltage. This is only true for low frequency AC signals. 
For high frequencies to ^ 'Tj^c^ the charge buildup can- 
not follow the changes in the AC voltage. This situation 
gives rise to an effective negative capacitance, or induc- 
tive behavior, since in this transport regime it appears 
as though the voltage on the plates lags the current, as 
shown in the parametric plot of the amplitude-phase di- 
agram in Fig. 7d. It should be noted that this effect 
arises from the presence of a quantum mechanical charge 
relaxation resistance, and is therefore quantum me- 
chanical in origin. Consequently, quantum mechanical 
devices exhibit inductive behavior at sufficiently high fre- 
quencies even though no geometric inductor is present in 
the circuit. 

It was further pointed out (271, that a second relevant 
time scale in high frequency transport of QDs is intro- 
duced by the carrier dwell time td of the QD, which 
is related to the line width function of the dot barrier 



7l = AH/td . Smaller 7^ serves the purpose of increas- 
ing the charge carrier dwell time inside the QD region 
and consequently lowering the transition frequency lutr 
at which quantum inductance appears. For sufficiently 
small 7l (i.e., large dwell times), the transport is always 
inductive, even at low frequencies; however, these regimes 
have not yet been realized experimentally due to the dif- 
ficulty that arises in fabricating strongly coupled topgate 
leads. Utilizing the above considerations, Wang et al., 
[131 define the quantum inductance as Lq = RqTd/l2. In 
essence, the finite line width forces a charge carrier to 
remain inside the dot for at least time to implying that 
the charge carriers can only follow changes in the voltage 
for frequencies w <C I/t^i. At higher frequencies this ef- 
fective trapping of the charges inside the QD for to gives 
rise to such large phase delays between the current and 
voltage that the capacitance appears negative. In our 
system, this dwell time effect is larger than the charge 
relaxation time, td > trc so that the inductive behav- 
ior is attributable to the trapping of charge carriers in 
the QD. 

B. Quantum Capacitance 

Unlike quantum inductance which is does not depend 
on graphene's particular density of states, quantum ca- 
pacitance, introduced by S. Luryi jsll, depends on the 
underlying band structure of the material. The quan- 
tum capacitance Cq describes the movement of the con- 
duction band as a function of the applied gate bias: 
Cq = ^ I ■ Unlike in most conventional semi- 
conductors, in graphene, the quantum capacitance is an 
important parameter since in the low bias regime, mono- 
layer graphene exhibits a linear gate tunable dispersion 
[22. 36, 37]. In the case of narrow constrictions, when 
the GNR width approaches the Fermi wavelength of the 
electrons, a bandgap opens (s^l and Cq (as well as the 
density of states) exhibit van Hove singularities, as shown 
in Fig. 2b [37|. Consequently the quantum capacitance 
of a graphene nanoribbon (per unit width) can also be 
tuned by varying the gate bias according to 

Cq {Vtg) = ^ I' 2 Q - ' (3) 

where 77 = Ep/ksT , a;„ = En/ksT is the energy of 
the n*'' subband, and O (77 — x„) is the Heaviside step 
function. The expression for the quantum capacitance of 
graphene stated in Eq. 3 is the same as the one obtained 
by other authors 1^37,], and it's derivation is outlined in 
Appendix II. We can combine the geometric capacitances 
of the gates and the quantum capacitance of the GNR to 
form the total capacitance 

11 1 

^Geometric 



4 




Figure 4: (a) real and (b) minus the imaginary components of the dynamic admittance of a 2.5 nm nanoribbon structure as 
shown in Figure la as a function of Vtg (Volt) and frequency, with Eo — /j. = 50mV, L = 200nm; (c) real (solid) and imaginary 
(dashed) components of the dynamic admittance for different values of the line width function 7^ = 0.1 meV (red), 7l — 0.3 
meV (blue), jl = 0.5 meV (green), using Vtg = 0.5V, using a 300 nm Si02 backgate and a 10 nm Hf02 top gate at T = 5 K. 



Since current graphene-based devices are fabricated using 
relatively thick backgate dielectrics (typically on the or- 
ders of hundreds of nanometers), CbGi consequently, has 
no contribution to C Geometric- We therefore only need to 
take into account the contribution from the topgate di- 
electric and can ignore the capacitive contribution from 
the backgate and the capacitive coupling to the contacts, 
which cannot be tuned by a gate bias. Although the ca- 
pacitive coupling to the backgate and contacts changes 
the value of the geometric capacitance, which in turn 
decreases the voltage range over which Cs can be modu- 
lated (a similar effect is achieved in Fig. 2a by changing 
the dielectric constant of the top gate), gate modulation 
of the total capacitance has been experimentally observed 
[sol . l40j . The total capacitance Cs is then determined 
by, Ce = Cq+Ctg ' '^^^ contribution of the quantum 
capacitance dominates in the limit of small quantum ca- 
pacitance Cq, i.e., in devices with thin gate dielectrics 
since it is in series with the geometric capacitance, and 
for capacitors in series, the smaller capacitance domi- 
nates p6l . [stI . m, |4l|. In typical GNR devices the oxide 
capacitance is on the order of C oxide ~ llbaF/iim^ 
and the trap capacitance, Ctr, which arises from the cou- 
pling to the leads is on the order of Ctr ~ IQfF/^m?. 
while the quantum capacitance can be tuned in the range 
5—35fF/iim? [3^. Figure 2a shows the total capacitance 
as a function of several typical dielectric materials. The 
quantum capacitance of a 2.5 nm GNR is shown in Fig 2b. 



While optical properties and carrier transport in GNR's 
are in general affected by the particular graphene edge 
state such as armchair, zigzag, or mixed edges (4^ . for 
simplicity we consider GNRs with pure armchair bound- 
aries [4^ . 



C. Fano Factor 

In order to compare our model with previous experi- 
mental work, we note that the real part of Eq. 1 is pro- 
portional to the frequency dependent shot noise [44, 4^ . 
Consequently, our model can be utilized to calculate the 
frequency dependent Fano factor of a graphene nanorib- 
bon. The dependence of the Fano factor of a GNR on 
gate and source-drain bias at w = 800 MHz is plotted 
in Figure 3. The device parameters, such as the di- 
electric thickness, GNR aspect ratio {W/L — 25, with 
L = 200nm) and operating frequency were chosen to 
match existing experimental results p3L Is^ l . The resul- 
tant dependence of the Fano factor is in excellent agree- 
ment with the measured values. Experimentally, the 
Fano factor peaks at a value J' — 0.34, while in our 
model it peaks a.t — 0.46. Furthermore, the change in 
slope of the Fano factor occurs at w 0.15 in both the 
experimental data and in our calculation. Although, in 
our calculation the Fano plot is symmetric about Vbg = 0, 
this is not the case in the experimental data. This dis- 
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crepancy is caused by the shift of the Dirac point in the 
presence of charged impurities [4^ |43|. The shght dis- 
crepancy between Fig. 3a and the measured value at 
zero gate bias and the theoretical prediction most likely 
stems from a finite density of states at the Dirac point, 
which is a feature measured in numerous experiments 
(isj . Although the above simulations model a relatively 
wide GNR, since the real part of the frequency depen- 
dent admittance, is proportional to the frequency depen- 
dent shot noise, the agreement between the experimental 
data and our model point to its validity. However, fur- 
ther experiments are needed for smaller GNRs at higher 
frequency to fully assess the theoretical predictions. 



III. RESULTS & DISCUSSION 



T = 0.01 meV 
T = 0.10 meV 
T = 0.20 meV 




0.3 0.4 
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Figure 5: Frequency dependent impedance Z (ixi) in the limit 
^ 0^ at various temperatures for the same device parame- 
ters as in Fig 4 (Vtg = 50 mV), and the same parameters as 
in Fig. 4. 



We now turn to the discussion of our main result, the 
dynamic admittance of the top gated GNR. The real 
and imaginary components of r(w) are plotted in Fig- 
ure 4. Figure 4b, shows the imaginary component of 
the admittance as a function of frequency and the top 
gate bias (the sign of the imaginary admittance has been 
inverted for clarity) where one can see that — 3r(a;) be- 
comes negative, corresponding to an inductive behavior 
of the device. The spikes in 3r (w) correspond to reso- 
nances between the electronic states inside the QD and 
the van Hove singularities of the GNR. Off-resonant cases 
are plotted if Fig. 4c. In general, the width of the GNR 
changes the energy spacing of the van Hove singularities. 
For a GNR with W=2.5nm, the spacing is ~1V (see Fig. 
2b), while for 5 nm GNR the energy spacing decreases 
to ~0.4V [37]. Consequently, in order to achieve the de- 
sired resonance source-drain voltages below the energetic 
spacing of the van Hove singularities should be selected, 
in our case Vsd — SOmy. Larger source drain volrages 
(or equivalently wider GNRs) would result in several van 
Hove peaks being inside the transport window. If a suf- 
ficiently small number of resonance levels (2 or 3) are 
present in the transport window, the sharp resonances in 
the conductance are no longer visible, however the gen- 
eral features of the complex admittance are still present. 
For a larger number of resonant levels, the transport be- 
comes diffusive, and the effect of quantum inductance is 
no longer observed. 

In the limit w ^> 0+ and large 7^, the real part of 
the frequency dependent impedance Z {uj) ~ (5RF)^^ — > 
hjle^ approaches the charge relaxation resistance Rq as 
shown in Figure 5. For small values of 7^, Z (uj) diverges 
since low values of the line width function correspond 
to large values of the carrier dwell times td, for which 
transport through the QD becomes blocked. 

For frequencies greater than wtj? where SF (uS) > 0, 
the transport through the GNR crosses over from being 
purely capacitive to RLC behavior and the admittance of 
the device can modeled by a classical series RLC circuit 
as 



^ ' 1 - uj'^LqC - iujCRq ■ 

The comparison of the classical RLC model to the quan- 
tum model is shown in Fig. 6. In the above model, 
the circuit inductance of the GNR is given by Lq — 
/i^/127re^7L and C — C^(0). In accordance with the 
model proposed by Wang et al. for a generic quantum 
dot [27] . the transition frequency ujtr in our device oc- 
curs below the resonance frequency wq. The transition 
frequency is determined by the carrier dwell time tjj (or 
equivalently jl) and ranges from 50 - 200 GHz for the 
values of 7^ shown in Fig. 4, which corresponds to in- 
ductances of 40 - 200 uH. Furthermore, the resonance 
frequency ujq of the resultant RLC circuit is given by 
ujQ — (LqC) . Using the above inductance values, for 
typical quantum dot sizes of 100 nm^ and 25 nm^ fabri- 
cated on 300 nm Si02 with a capacitance of 11.5 nF / 
cm^ and 30 /iF/cm^ quantum capacitance of graphene, 
we predict resonant frequencies of this oscillator on the 
order of 100 - 900 GHz. This gives rise to Q-factors 

Q= ft = Ti \[h ™ ^^'^^^ of Q = 1 X 10^ - 8 X 10^ 
(Aw is the width of the resonance) . This is a remarkable 
number for an all electronic circuit. To compare this re- 
sult to conventional 2DEG structures, such as AlGaAs 
hetero junctions, we note that the dwell time may be ap- 
proximated as Tu ^ 4:L/vf [231, where L is the device 
length and vp is the Fermi velocity. Since vp is AlGaAs 
2DEGs is typically on the order of 3 x lO^m/s, the cor- 
responding Quantum inductance and Q factor would be 
an order of magnitude lower. It should be noted that the 
above results were obtained assuming prestine nanorib- 
bons, and effects that would lead to dephasing or the 
degradation of the ballistic mean free path, such as the 
morphology [49| of the device and potential fluctuations 
have not been considered and are beyond the scope of 
this study. In general, spatial potential fluctuations in 
graphene are on the order of tens of meV [8], while the 
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Figure 6: A Comparison of (a) real and (b) imaginary admit- 
tance components of the classical RLC model (solid black line) 
with the results of our simulation (dotted red line) 7l =0.1 
meV and 7l = 0.3 meV for a 200 x 200 nm quantum dot 
(which is easily fabricated using modern technology), and the 
same parameters as in Fig. 4. 



pie, if a high frequency AC signal is applied to the source- 
drain contacts, a subsequent sweep of the top gate voltage 
would drive the oscillator in and out of resonance, allow- 
ing for an all electronic measurement of this frequency 
via the device gain. Consequently, a measurement of the 
on and off resonant gain in the device yields an all elec- 
tronic measurement of the frequency of the signal. Since 
Wo is on the order of 100+ GHz, this GNR-RLC circuit 
could be used to measure frequencies which are currently 
unattainable by standard instrumentation. This, coupled 
with graphene's intrinsically ultra-high carrier mobilities, 
render it an ideal material for all-electronic THz devices. 



IV. CONCLUSION 

In conclusion, we have applied the Green-Kubo for- 
malism to model high frequency transport through dual 
gated graphene nanoribbons. We showed that above a 
sufficiently high frequency WTi? determined by the dwell 
time of charge carriers in the gate-defined quantum dot, 
the behavior of the GNR device is analogous to a classical 
RLC oscillator with a very high Q-factor. The inductive 
behavior arises from the negative capacitance of the QD 
which occurs when charge carriers become trapped inside 
the dot for times r^j and cannot follow the driving volt- 
age. This leads to a phase lag greater than tt and thus 
results in a negative capacitance. Coupling the inductive 
behavior of a quantum dot to the gate tunable quantum 
capacitance in graphene, gives rise to an in-situ tunable 
ultra high frequency oscillator and filter thereby extend- 
ing the reach of high frequency electronics into the THz 
regime. 
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Appendix I 

Here we derive the complex admittance for a QD with 
a single resonant level. We emphasize that these results 
are general for any QD and are not specific to graphene. 
Following the procedure in [s^l, the expression in Eq. 1 
for the admittance may be expressed in terms of trans- 
mission T (E) and reflection R (E) amplitudes as 



2e2 



h 2t. 



dEi 



dE^Ii^^t-Ii^x 
uj + E-i ~ E2 + ie 



7 



3? {T (El) T* {E2) + l-R{Ei)R* {E2)} , (4) 

where / is the Fermi function. The transmission and 
reflection amphtudes through the barrier are given by 



T{E) = 



E-E0 + i-fL/2 



R{E) 



E — En 



E-Eo + i-/l/2' 



where 7^ is a hnewidth function representing the cou- 
phng of the lead to the QD and characterize the energy 
width of the resonance |30|. The hnewidth function of 
the lead is related to the carrier dwell time, tjj, inside 
the QD according to 

Ah 

IL = • 

TD 

Substituting Eq. 5 into Eq. 4 and carrying out the inte- 
grations, we obtain the following result for the dynamic 
impedance of the QD T (w) = (w) + (w): 
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Figure 7: (a) Real and (b) Imaginary components of the quan- 
tum dot admittance (which does not include the contribution 
of the GNR) Fo (1^); (c) several cuts through (a) at = 2 
(dotted black line), C, — 0.2 (dashed blue line) and (b) at 
C ~ 2 (dot dashed green line) and ^ = 0.2 (solid red line); 
SiF_D {(jj) > corresponds to an inductive behavior; (d) para- 
metric plot of real and imaginary components as a function 
of to for the same values of ^ as in (c); 



G(^) = ^ / dET{E)M{E)-l^, 



where /g is the Fermi function. In the ballistic limit 
T (E) — s> 1. Furthermore, the mode function M (E) is 
given by 



In 



(C7l/2 + M' + (7l/2)' 



(C7L/2 - hLof + (7l/2)' 



(C7L/2)' + ilL/2y 



(6) 

The real and imaginary components of F/j (uj) are plotted 
in Figure 7 as a function of the dimensionless parameter 
C (see main text for discussion). 



M 



where Vg is the group velocity given hy Vg (^k 
\-§§-- Using this, we can express the mode function as 



M{E) 



WL 

47r2 



TT 00 



de kdkS E-E k 



hvpkx 



which becomes. 



Appendix II 



In this section we provide a brief derivation of the 



quantum capacitance, Cq = 



dE\E=EF' 



To do so we 



require an energy dependent expression for the density of 
states. This is computed in the Landauer-Biittiker for- 
malism. In the low bias regime the ballistic conductance 
(^(s) jg given by 



M{E) 



2W \E\ 



where the factor 2 is introduced to account for spin 
degeneracy. Inserting the above expression for M{E) 
into G^'^\ and making use of the fact that — 



g(E-EF)/fc£ 



[1+e 



{E-Ep)/kj: 



-2-, we obtain 



8 



2e^ 2W 
h hnvp 



dEE- 



^(-E~EF)/kBT 



(l + e(~-E--E^)/'=«^)' 



These results can now be used to calculate the quan- 
tum capacitance of a graphene strip defined as Cq = 
I e=Ef ' Following the same procedure as above, the 
carrier densities are given by 



dEE- 



^{E~EF)/kBT 



{IJ^ e'^E-EF)/kBTy 



Integrating the above expression by parts, and doing 
some algebra, we obtain 



2P 

''-P= -J^r (2) { J-i {-f^p) + ^1 {rip)} , 

TT n U rp 



which is in agreement with literature [37]. Using this 
expression, and differentiating it at the Fermi energy, we 
obtain 



G(^) = ^^|M:r (1) {^0 {-vf) + -Fo , 
n nirvF 

where T (x) is the Gamma function, J-s (rj) is the 



Fermi-Dirac Integral given by J-^ {r]) 



s+l) I 



r(s+l) J ■ 




and rjp — -r^- Taking the low temperature limit we 



have 



-dfo 
dE 



S {E — Ep), in which case the expression 



for the ballistic conductivity reduces to 



(e) = (2) { j-o i-VF) + To (vf)} , 



which in the low temperature limit, can be shown to 



be 



Cr 



2e^kBT 



In 



2 1 + cosh 



qVTG 



G 



(B) 
T=Q 



2e^ 2W\Ef\ 
h ttHvf 



This result can be further generalized to a case of a 
realistic GNR of finite width, yielding Eq. 3 in the text. 
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